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Abstract 

In this letter, we demonstrate a strong dependence of the electrostatic deformation of doubly- 
clamped single-walled carbon nanotubes on both the field strength and the tube length, using 
molecular simulations. Metallic nanotubes are found to be more sensitive to an electric field than 
semiconducting ones of the same size. For a given electric field, the induced deformation increases 
with tube length but decreases with tube radius. Furthermore, it is found that nanotubes can be 
more efficiently bent in a center-oriented transverse electric field. 

PACS numbers: 85.35.Kt, 32.10.Dk, 31.15.Qg 



1 



Carbon nanotubes (CNTs) can be used in nanoelectromechanical systems (NEMS) for 
uses in sensing, actuation, vibration, and laboratory-on-a-chip applications pQ . Particularly, 
doubly-clamped (suspended) structures of CNTs have been reported to be used as key 
components in a number of nanodevices [21 EJ IH EJ El [7] . In most of these recent NEMS, 
nanotubes are usually suspended between two electrodes with an applied voltage. However, 
as to be demonstrated in this paper, even without applying electric current and excess 
charges, doubly-clamped CNTs can still be significantly bent in a static transversal electric 
field, as a result of electric polarization. 

Cantilevered CNTs have long been known to deform in an electric field [SI El HO] • The 
mechanism of this deformation relies on the fact that CNTs can be bent by a bending 
moment induced by interactions between the electric field and the molecular dipoles due to 
electric polarization. This property has been exploited in the design of nanorelays [HI [12] . 
Compared with cantilevered CNTs, doubly-clamped ones can have many advantages in 
electronic devices [13] (e.g. they can be integrated in nanodevices with a well-defined spatial 
structure and resonance frequency). This letter investigates the deflection of neutral single- 
walled CNTs (SWCNTs) by an electrostatic field. Two semiclassical theories have been 
combined for characterizing chemical potential and electrostatic interactions in both metallic 
and semiconducting CNTs. Considering the general correlation between the conductivity 
and polarizabilities of CNTs, metallic CNTs are supposed to be more sensitive to electric 
fields than semiconducting ones are, and hence can be expected to play a more important 
role in NEMS. 

When an electric field is applied to a metallic CNT, an induced dipole can be created at 
each atomic site with a quantity of free charge, by shifting the electrons and the nuclei. In 
our calculation, each atom is therefore modeled as an interacting polarizable point with a 
free electric charge, while the chemical bonds are described by using a many-body potential 
function. Motionless equilibrium positions of carbon atoms in an electric field are computed 
by minimizing the total potential energy of systems U tot , which is the sum of two terms: 
jjtot _ jjeiec _|_ jjp^ w ;h ere jjv i s th e interatomic potential due to the C-C chemical bonds 
in absence of an external field, including the long-range interactions. JJ elec stands for the 
electrostatic energy from the interaction between charges, dipoles and external fields. JJ elec 
is calculated using a Gaussian-regularized charge-dipole JTU] and a dipole-only [H] model 
for metallic and semiconducting CNTs, respectively. U p is computed using an adaptive 



2 



interatomic reactive empirical bond order (AIREBO) potential [15J, which is an extension 
of the REBO potential [TB] which includes a Lennard- Jones potential. Compared to first- 
principle and semi-empirical methods, an important feature of our combined models is their 
ability to deal with large systems (up to 5000 atoms). This is particularly important for 
studying geometric effects on the electrostatic bending of CNTs, since the axial periodic 
condition can hardly be applied in this issue. 

The strength of electric field E used in this work is in the order of V/nm, because the 
tubes used in our calculation are too short (18 ^ 27 nm) to be bent in a weak electric field. 
As will be shown, the required field strength decreases with increasing tube length L for 
a given deformation. Considering that the CNTs used in experiments are usually of /im- 
length, field emission effects p2] and conductance switching due to strong transversal fields 
[18] are both neglected. These field strengths are not large enough to cause a significant 
change in the chemical bond strength [19]. Moreover, it can expected that an armchair CNT 
can exhibit more deformation by an electric field than a metallic zigzag tube of the same 
size, since previous studies concluded that bending deformation would have a negligible 
effect on the electronic transport properties of armchair CNTs [201 El] , while the quantum 
conductance of zigzag CNTs can significantly decrease under large bending deformations 
[22]. In this work, the tubes are assumed to be suspended between two electrically insulated 
supports. 

The mechanism of electrostatic bending of a SWCNT is depicted in Fig. [T] (a). When a 
CNT (initially electrically neutral) is submitted to a transversal electric field E ext , the field 
tends to shift negative and positive charges in opposite directions (in this figure we can see 
that positive charges move to the top and negative ones move to the bottom of the tube), 
two induced molecular dipoles (pi and P2) are hence created in the tube. The tube is then 
bent by two bending moments (tl and r 2 ), which is induced by the interaction between the 
molecular dipoles and E ext (r = p m <g> E ext ). 

Molecular dynamics (MD) simulations were performed using the AIREBO potential, for 
generating arbitrary initial configurations of CNTs at room temperature. The energy of these 
MD-generated configurations is minimized in order to calculate their equilibrium positions 
in an electric field. MD simulations showed that CNTs slightly oscillate around their initial 
positions due to thermal fluctuations, which have been observed in experiments [2"B"j 12%]. We 
found that the orientation of electrostatic deformation of CNTs strongly depends on their 
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FIG. 1: (Color online) Schematics of electrostatic deformation of a SWCNT. (a) Two arrows p\ 
and P2 stand for two induced molecular dipoles. E ext is represented by the arrows with dashed 
lines, n and t<i stand for the induced bending moments acting on the tube. The color scale of 
atoms is proportional to the density of induced charges. The vectors stand for the induced atomic 
dipoles. (b) Two possible deflection directions in bending (upward or downward). 

initial position at the moment when E ext is applied, e.g., as shown in Fig. [I] (b), if E ext is 
applied when the tube is at position 1, its equilibrium state in electric field will be at position 
3; conversely if the tube is initially in position 2, it will be bent to position 4. These results 
also imply that no change of the deformation direction will take place if the field direction 
is totally reversed, due to the system symmetry. For measuring the amplitude of deflection, 
we define u as the displacement of the center of the tube middle from its initial position. 
For a small electric field, the bending of the CNT is elastically reversible, i.e., the tube will 
come back and oscillate around its initial position once E ext is removed [HI |25] . In our MD 
simulations, the average frequency of thermal oscillation (without electric field) of a (5, 5) 
tube (length L « 20.0 nm) is found to be 60 ± 10 GHz for the several first harmonics. This 
value is comparable to oscillation frequency of doubly-clamped CNTs recently reported in 
Ref. [26]. 

We plot in Fig. [2] computed values of u versus external field strength E for both a (4, 
4) and a (5, 3) CNTs. It can be seen that the metallic tube is clearly more sensitive to 
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FIG. 2: (Color online) Electrostatic deformation it versus field strengths E = |.E e:c *| for two 
SWCNTs: a metallic (4, 4) (L ~ 24.5 nm) and a semiconducting (5, 3) with almost the same 
length and radius. The tube (4, 4) begins to break down when E > 11.0 V/nm, then failure occurs 
very soon thereafter {u max > 5.1 nm). The inset shows the position of fracture and the distribution 
of von-Mises stress at the bottom (before the fracture occurs). 

E ext than the semiconducting one due to their different polarizabilities [27J. When u ~ 5.5 
nm, it is observed that the (4, 4) tube begins to fracture near the fixed end, in which 
the maximal tensile strain is localized (inset in Fig. [2j. We can also see that there is no 
deformation when E remains small, since in such cases the moment of electric force are not 
large enough to exceed the mechanical resistance barrier of the tube. These results also 
imply that semiconducting tubes can sustain higher field strengths. 

The CNTs used in experiments are usually microns in length. However, their field-induced 
deformations can hardly be addressed by the calculations using atomic models in a direct 
way. Hence, it is necessary to address the influence of the tube length L. In Fig. [3] (a), it 
can be seen that the relative displacement u/L increases with L, i.e., a weaker electric field 
is required for a given u in a longer tube. This increased displacement depends on the field 
strength. An extrapolation was done from the data E versus L for tube lengths up to 1 /im, 



5 




J 1 1 1 1 1 1 

20 25 30 35 40 45 50 

L (nm) 

10 -] 1 



I - 

> 

LU 

0.1 - 

15 100 1000 

L (nm) 

FIG. 3: (a) u/L versus the tube length L for (4, 4) SWCNTs, in 4 different transverse E ext (which 
is perpendicular to the initial axis of the tubes), (b) E versus L on a logarithmic scale for two given 
u. The symbols present the calculated points and the lines stand for the extrapolation curves. 

in order to make experimental verification possible (Fig. [3] (b)). It was found that ln(E) 
decreases almost linearly with ln(L) for two given u. This is related to the fact that the 
electric capacitance per unit length of a metallic cylinder is roughly proportional to ln(L). 
Thus, we can conclude that, for the //m- length tubes used in experiments, the required field 
strength is on the order of V/ /im. 

To explore further geometry effects, we plot u versus E for the tubes of different radii 
R in Fig. |4| It is found that small tubes are more sensitive to E ext than large ones. This 
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FIG. 4: u versus E for 6 individual armchair SWCNTs (L ~ 24.5 nm) with different radii. The 
direction of E ext is perpendicular to the initial axis of tubes. Inset shows the force balance between 
an electric moment and internal stresses in a half of a deformed tube. 

looks abnormal since it is known that the polarizability of a CNT increases with its radius. 
However, from a mechanical point of view, the tube becomes more difficult to be bent due to 
the increase of its area moment of inertia (J = 0.5m(Rj nt + R^ ut ) for a thick-walled cylindrical 
tube, where m is the mass). From our results it is clear that the later effect plays a more 
important role. The ratio u/E is roughly the same for all of these tubes in case of large 
deformation (u > 1 nm). 

In Figs. H and g] we can see a threshold field of a few V/nm for each u versus E curve 
before one obtains a nonzero deformation. The existence of this threshold field relies on the 
fact that the electric bending moments r remains weaker than the tube resistance when E 
is small. With increasing E, one reaches the limit of the threshold field until a value of u 
allows a moment balance in the nanotube (inset in Fig. [4]). It is found that the threshold 
field increases with R but decreases with L. 

In this study, the applied external field is assumed to be uniform as a common theo- 
retical simplification. However, in experimental situations, the external fields are usually 
not homogeneous and strongly depend on the experimental set-up. By changing the field 
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FIG. 5: (Color online) (a): Schematic (in y-z plane) of the deformation of a CNT subjected to an 
y \-like electric field E ext (dashed lines). The field angle 6 is defined as the angle between the field 
lines and the axis z. The field vector E = (E x , E y , E z ) = (0, Esm(9), +E cos(8)) when z < 
and E = (0, Esm(0), -Ecob(6)) when z > 0. (b): u versus E for a (4, 4) SWCNT (L m 24.5 
nm) in E ext with different 6. The simulation data are represented by the symbols with eye-guide 
lines. 

direction, we found that CNTs can be more efficiently bent in a center-oriented (f \-like) 
electric field (see Fig. [5] (a)). Results of the deformation of a SWCNT by this field is plotted 
in Fig|5] (b). We can see that values of u roughly follow a linear relationship with E when 
u > 1 nm. It was also found that, when 9 = 7r/4, the tube can be most efficiently bent and 
its u versus E curve has straightest shape. 

In summary, deformation of doubly-clamped SWCNTs in an electrostatic field has been 
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simulated using a charge- dip ole polarization model combined with an empirical potential. 
The interplay between the mechanical resistance and electric polarization of CNTs was 
investigated for both metallic and semiconducting tubes. Metallic CNTs are found to be 
more sensitive to an electric field than semiconducting ones. This study reveals that the 
field-induced deflection increases with tube length while it decreases with tube radius. It 
was also found that CNTs can be more efficiently bent in a center-oriented electric field. 
From a theoretical point of view, graphene nanoribbons and other metallic nanowires/tubes 
should have similar properties, depending on their dielectric constants. This electrostatic 
deformation of doubly-clamped nanostructures can be expected to open a path to designs 
of novel nanodevices. 
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